SPHERICAL MEANS IN ODD DIMENSIONS 
AND EPD EQUATIONS 
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^ ' Abstract. The paper contains a simple proof of the Finch-Patch- 

Q , Rakesh inversion formula for the spherical mean Radon transform 

1/'^ ' in odd dimensions. This transform arises in thermoacoustic to- 

mography. Applications are given to the Cauchy problem for the 
Euler-Poisson-Darboux equation with initial data on the cylindri- 
cal surface. The argument relies on the idea of analytic continua- 
tion and known properties of Erdelyi-Kober fractional integrals. 
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c^'. 1. Introduction 

We consider the spherical mean Radon transform 

., (1.1) (M/)(^,t) = ^/ f{e-ta)da, 

g: ^g5"-\ tGM+ = (0,oo), 

where / is a smooth function supported inside the unit ball B = 

{x G M" : |x| < 1}, S""'^ is the unit sphere in M" with the area 

cr„-i = 27r"/^/r(ri/2), and da denotes integration against the usual 

O ■ Lebesgue measure on S"""^. The inversion problem for this transform 

has attracted considerable attention in the last decade in view of new 
developments in thermoacoustic tomography; see |AKKt AKQ, IFPRt 



IFHRl IKKt IKut IPSll IPS2J and references therein. Explicit inversion 
5t , formulas for Mf in the "closed form" are of particular interest. For 

n odd, such formulas were obtained by Finch, Patch, and Rakesh in 
|FPRj . The corresponding formulas for n even were obtained by Finch, 
Haltmeier, and Rakesh in |FHRj . An interesting inversion formula of 
different type, that covers both odd and even cases, was suggested by 
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2 B. RUBIN 

Kunyansky |Ku] : see also a survey paper [KKj . where diverse inversion 
algorithms and related mathematical problems are discussed. 

In spite of their elegance and ingenuity, most of explicit inversion for- 
mulas for Mf are still mysterious, and the basic ideas behind them are 
not completely understood. We observe, for example, that the deriva- 
tion for n = 3 in [FPR] relies on implementation of delta functions 
and is not completely rigorous. The dimensions n = 5, 7, . . . have been 
treated using a pretty complicated reduction to the 3-dimensional case. 
The resulting inversion formula can be written in the form 



(1.2) /(x)=c„a/ [D""3r-Ve] 

where ipe{t) = {Mf){e,t), 



t=\x-e\ 



■^19x1' """"4^2-11(71/2)' 2tdt' 

cf. |FPRt Theorem 3]. Other formulas in that theorem can be obtained 
in the framework of the same method; see Remark 13. 1[ 

In the present article (Section 3) we suggest a simple and rigorous 
proof of 01. 2p that handles all odd n simultaneously, so that no reduc- 
tion to ra = 3 is needed. The idea is to treat the spherical mean as 
a member of a certain analytic family of operators associated to the 
Euler-Poisson-Darboux equation 

n + 2a — 1 

(1.3) UaU = A^U-Utt Ut = 

and invoke known facts from fractional calculus [SKMj about Erdelyi- 
Kober operators. Section 4 deals with applications. Here we recon- 
struct the solution u{x,t) of the equation D^u = X^u, A > 0, from the 
Cauchy data on the cylinder 5*""^ x M+. Section 5 contains comments 
and open questions. 

No claim about the originality of the results presented in this article 
is made, but it is felt that the elementary use of operators of fractional 
integration to obtain them might appeal to the applied mathematician. 
If the reader is not interested in applications in Section 4, he may skip 
subsections 2.2 and 2.4. These are not needed for the basic Section 3. 

Acknowledgements. I am grateful to professors Mark Agranovsky, 
Peter Kuchment, and Larry Zalcman for very pleasant discussions and 
hospitality during my short visits to Bar Ilan and Texas A&M Univer- 
sities. 
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2. Preliminaries 

2.1. Erdelyi-Kober fractional integrals. We remind some known 
facts; see, e.g., jSKMl Sec. 18.1]. For i?ea > and r/ > -1/2, the 
Erdelyi-Kober fractional integral of a function yj on R+ is defined by 

(2.1) (/»(t) = / {f - r'r-'r'^^'^ir) dr, t > 0. 

r(a) Jq 

For our further needs, it suffices to assume that ip is infinitely smooth 
and supported away from the origin. Then I°'ip extends as an entire 
function of a and r], so that 

(2.2) Jjy. = v., 

(2.3) (/»(t) = t-'("+'')Z}"^t'(°+"^+'')(/°'""<^)(t), ^ = ^ T ' 

(2.4) {l;;""f)it) = t-^^'^'"'^D"'f'^if{t), 
where m is a nonnegative integer. The property 

allows us to write (12.31) and (12. 4p in a different equivalent form. The 
composition formula and the inverse operator are as follows: 

(2.6) j,v> = ^r^"^, (^.")"V = /,-+„^- 

2.2. The generalized Erdelyi-Kober fractional integrals. Let J^ 

and Ii, be the Bessel function and the modified Bessel function of the 
first kind, respectively [EJ. The generalized Erdelyi-Kober operators 
are defined by 

(2.7) j,",y.(t) = t~2("+^) j,"t2v(t), i:;Mt) = t-'("+^)/rt'v(t), 

where A > 0, 



(2.8) J^ifiit) = 2°A^-" / it^-r^Y''-'^^^Ja-i{\Vt^^^)ip{r)rdr, 

Jo 

(2.9) JXt) = 2''X^-'' [ {t^-r^Y''-^^^^I^^i{\VW^^)ip{r)rdr; 



JO 

see |L11IL2] . |SKMl Sec. 37.2]. As above, we assume ip to be infinitely 
smooth and supported away from the origin. Integrals (12. 8p and (12. 9p 
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are absolutely convergent if Re a > and admit analytic continuation 
to all complex a by the formulas 

(2.10) J> = D"''J^+"'^ = J"+"'D"'^, 

(2.11) /> = D™/°+"^(^ = /°+"^D'^y?, 

m G N. These follow from the well-known relation 
1 d^"^ 



r CLT ^ 

(similarly for Q. Clearly, J^q = /"q = I^. liRea>0 and Rei3>0, 
then, changing the order of integration and using the formula 2.15.15(1) 
from [PBMj . we get 

(2.12) Ca'^.> = CV 

The latter extends by analyticity to all complex a and /3 and yields the 
inversion formula 

(2.13) iJ-,r'f = i;;:^j. 

By (12.111) and (12. 7p . this can also be written as 

(2.14) (J^"a)"V = t-2'?Z}™/™-"t2{^+a)^ _ ^-2r,jjrn^2iv+m)jrn^aj_ 

2.3. The Euler-Poisson-Darboux equation. Consider the Cauchy 
problem for the Euler-Poisson-Darboux equation (II. 3p : 

(2.15) n^u = 0, u{x, 0) = /(x), ut{x, 0) = 0, 



where / belongs to the Schwartz space iS(]R"'); see |B1] for details. If 
a > (1 — n)/2, then (12.151) has a unique solution u{x,t) = {M^f){x) 
where the operator M" is defined in the Fourier terms by 

[M^fny)=m^it\y\)fiy), 

m^{p) = r(a + n/2) (p/2)i-°-"/2 J^^,^^_,{p). 

The operator M" extends meromorphically to all complex a with the 
poles —n/2, —n/2 — 1, ... . For Re a > 0, it is an integral operator 

r{a + n/2) , ,^ ,^,,2,,_i 



(2.16) (Mf /)(a;) = \,,^/ ' / (1 - IvlT-'fix - ty) dy. 
In the case a = 0, M"/ is the spherical mean 

(2.17) (M°/)(x) = ^ / f{x- ta) da; 
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cf. (11 -ip . Passing to polar coordinates, one can obviously represent 
M"/ as an Erdelyi-Kober integral of the spherical mean 

(2.18) (Mf/)(x) = ^^^±^(/>.)(t), ^.(t) = (M°/)(x), 

with 1] = n/2 — 1. 

2.4. The generalized Euler-Poisson-Darboux equation. Consider 
the more general Cauchy problem 

(2.19) n„M = AV u{x, 0) = fix), ut{x, 0) = 0, 

where / is a Schwartz function and A > 0. If a > (1 — n)/2, then (12.191) 
has a unique solution u{x,t) = (M";^/)(x), where the operator M^-^ is 
defined as analytic continuation of the integral 

(2.20) (M^/)(.)^ ^^/^)"";;r^^^^ ^— 



fix - y)it' - |2/n("-^)/2 Ja-ii\V^^W) dy; 

\y\<i 
see |B2] for details. As above, 

(2.21) (M^,/)(a;) = ^^^+^y\ j,>.)(t), ^.(t) = (M°/)(x), 

7] = n/2 — 1, where J^^ is the generalized Erdelyi-Kober operator (12. 7p . 
2.5. More preparations. We restrict M"f to x = 9eS'^^^ and set 

(2.22) {N''f){e,t) = e+^''-^{M^f){e). 

In particular, for i?ea > 0, owing to (12.161) . we have 

(2.23) (iV"/)(g,t) = ^^; + ff / f{y){e-\y-e\Y^-^dy 



where (...)"~^ has a standard meaning, namely, (a — 6)"^^ = (a — 6)°^^ 
if a > 6 and otherwise. Given a function F on the cylinder S*"^^ x M+, 
we denote 

(2.24) {PF){x) = ^—f F{9,\x-9\)d9, xeW, 

which is a modification of the back-projection (or dual) operator; cf. 
[HI |N], where these notions are used for the classical Radon transform. 
We also invoke Riemann-Liouville integrals [SKM] 

(2.25) {r^u){s) = -^ f {s-tr-'u{t)dt, ueC^[-lM 

r(«) J-i 
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The integral fl2.25p is absolutely convergent when Rea > and extends 
by analyticity to all a G C, so that 

(2.26) {IZ^'u){s) = {d/ds)'^u{s), m = 0,l,2,.... 

Lemma 2.1. Let B = {x & M."- : \x\ < 1}. For Rea > and any 
integrable function f supported in B we have 

(2.27) (PiV-/)(x) =cj ^^y\ ^ {I\u){h) dy, x e B, 

Jb \x-y\^ 

where 

_ Tin/2)Tia + n/2) ,J-\^-\y\\ ,(,) = (i _ ,.)(.-3)/.. 



-« 2i-°7r(«+i)/2r((n-l)/2)' 2\x-y 

Proof. By (12.231) and (12.241) . changing the order of integration, we get 

{PN-f)ix) = ^^" trff , / de [ fiy)i\x-e\'-\y-e\%-Uy 

cr„_i 7r"/2 T{a) Js„-i J^ 
T{a + n/2) 



where 



f{y)kaix,y)dy, 



ka{x,y) = [ {\x\^-\y\^-2e-{x-y))l-Ue 



(^n-2 



(/l-t)°-l(l-t')("~3)/2^^_ 



{2\x-y\y- 
This gives the result. D 

3. Inversion of the Spherical Mean for n Odd 

Let C°°{B) be the space of C°°-functions on M" supported in B] 
f e C^{B). By ^2?m . (I2A8D . and ([MD, analytic continuation of A^"/ 
at a = 3 — 72 has the form 

{N'-'m.t) = ^^Y(^^^"''^""V.W, Mt) = iMf)ie,t). 

Iin = 2k + 3, k = 0,l,..., then JKm yields 

(/!Y"n)(/i) = {d/dh)''\i - ef = {-ifk\ = (-i)("-=^)/2r^^ _ 2). 

Hence, analytic continuation of (12.271) at a = 3 — n is 

(3.1) P[D"-3r-V.](x)=c(/V)(x), c = 2(-l)("~3)/2r2(n/2)/vr, 
where x & B and 

. N . •? r^. ^ r(n/2-l) f f(y)dy 

(3.2) {I f){x) = ^ ' ' ' ■'^^' ^ 



47r"/2 Jj^ \x - y|"-2 
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is the Riesz potential of order 2 [SKMj . The latter can be inverted by 
the Laplacian, and simple calculations yield 



(3.3) /(a:)=c„A / [D^-V^Ve] 



(-l)in-l)/2 



t=M| ' " 4W2-ir(n/2)' 



Remark 3.1. Formula (13.31) coincides (up to notation) with the third 
formula in [FPRt Theorem 3]. Other formulas in that theorem can be 
similarly obtained from (12.271) if the latter is applied to A/ instead 
of /. Here we take into account that PAf = — /, because supp/ is 
separated from the boundary of B. 

4. An inverse problem for the EPD equation 

Let a > (1 — n)/2, A > 0. Suppose we know the trace u{9, t) of the 
solution of the Cauchy problem 

(4.1) n„M = X\, u{x, 0) = f{x), ut{x, 0) = 0, 

on the cylinder {{0,t) : 6 G S*""^, t G M+} and want to reconstruct the 
initial function / in the space C°°{B). This can be easily done using 
(13.11) and the Erdelyi-Kober operators. Indeed, by (12.211) . 

where (/^^(t) = (M°/)(e), r/ = n/2-l. Then by fl2A3ll . for ue{t) = u{9,t) 
we have 

^ _ T{n/2) _ T{n/2) 

Now, ([3IID yields 
nl+n/2) Pi^r't-%:^,,ue]=cl'f, c = 2(-l)(-^)/^r^(n/2)A, 



and therefore. 



(-l)in-l)/2 

V / A or ni— 3j.n— 2 r— o 



The operator I~^^x ^^"^ t)e replaced by any expression from (I2.14p . 

It remains to note that once / is known, the solution u{x,t) of the 
equation DaU = X^u can be reconstructed from the trace u{0, t) by the 
formula u{x,t) = {M^J){x); see ^TM . 
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5. Comments 

1. It is a challenging open problem to appropriately adjust our 
method to the case when n is even and give alternative proof of the 
corresponding inversion formulas from [FHRj and |Kuj . 

2. Formula (13.11) provokes the following 

Conjecture. Let n > 3 be odd. A function ip0{t) = ip{6,t) belongs to 
the range of the operator f —>■ {Mf){6,t), f G C°°{B), if and only if 
P[L)[^-^r-Ve] belongs to the range P[C°^{B)] of the potential / TOP . 

The "only if part follows immediately from (13.11) . The "if part 
requires studying injectivity of the back-projection operator P, which 
is of independent interest; cf. [Rj, where injectivity and inversion of 
the dual Radon transform is studied in the general context of affine 
Grassmann manifolds. Various descriptions of the range of the spher- 
ical mean transform and many related results can be found in |AKQ 

EH EH- 

3. It is worth noting that for n = 3, the inversion formula for Mf 
becomes elementary if / is a radial function, i.e. /(x) = /o(|x|). 

Lemma 5.1. /// G LL(K"), f{x) = fo{\x\), then {Mf){e,t) = Fo{t), 
where 

(5.1) Foit) = ,,,^,, ^ \ ,,, / /o(r) [a{r,t)r-'rdr 



7ri/2r((n-l)/2) 



|i-t| 



a{r,t) = [r^ — (1 — i)^]^''^[(l+t)^ — r^]-'^/^/4 being the area of the triangle 
with sides l,t,r. 

Proof. 

{Mf){d,t) = ^ I fo{\e-ta\)da 

= ^^^ f /,(Vl+t2_2ts) (1 - s2)("-3)/2 ds 
CTn-l J -I 

and (O) follows. D 

Corollary 5.2. If n = 3 and f is supported in B, then ( 15. 1\) yields 

1 '-' 



(5.2) (a) Fo(t) = -/ /o(r)rrfr if < t < 1, 

and 



'^t-n-t 



(5.3) 



ih) F,{t) = ^ [ foir)rdr if 1 < t < 2. 
^t^ Jt-i 
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In the case (a), 
(5.4) 



In the case (b), 



(5.5) 



Mr) 



fo{r) 



imt)) 



t=l-r 



(tFoit)) 



t=l+r 



Remark 5.3. We note that in (15. 4p and (15. 5p it is not necessary to 
know -Fo(^) for all t G (0,2) as in (13. 3p . It suffices to know it only 
for t G (0,1) or t G (1,2). It would be interesting to obtain similar 
inversion formulas for {Mf){6,t) in the general case, when [Mf){6,t) 
is known only for {e,t) G 5""^ x (0, 1) or {9,t) G 5"^"^ x (1,2), as in 
the radial case. 
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